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PREFACE 


The  work  described  in  this  report  was  authorized  by  Headquarters, 

US  Army  Corps  of  Engineers  (HQUSACE) ,  as  part  of  the  Coastal  Problem  Area  of 
the  Repair,  Evaluation,  Maintenance,  and  Rehabilitation  (REMR)  Research 
Program.  The  work  was  performed  under  Work  Unit  32415,  "Experimental  Testing 
of  Methods  for  Reducing  Wave  Runup  and  Overtopping  on  Structures,"  for  which 
Mr.  John  P.  Ahrens  was  Principal  Investigator.  Mr.  John  H.  Lockhart,  Jr. 
(CECW-EH) ,  was  the  REMR  Technical  Monitor  for  this  work. 

Mr.  William  N.  Rushing  (CERD-C)  was  the  REMR  Coordinator  at  the  Direc¬ 
torate  of  Research  and  Development,  HQUSACE;  Mr.  James  E.  Crews  (CECW-0)  and 
Dr.  Tony  C.  Liu  (CECW-EG)  served  as  the  REMR  Overview  Committee; 

Mr.  William  F.  McCleese,  US  Army  Engineer  Waterways  Experiment  Station  (WES), 
was  the  REMR  Program  Manager.  Mr.  D.  D.  Davidson,  Wave  Dynamics  Divi¬ 
sion  (WDD) ,  Coastal  Engineering  Research  Center  (CERC) ,  WES,  was  the  Problem 
Area  Leader. 

The  work  was  performed  at  WES ,  and  this  report  was  prepared  by 
Mr.  Donald  L.  Ward  and  Mr.  Ahrens,  WDD,  CERC,  under  the  general  supervision  of 
Dr.  James  R.  Houston,  Director,  CERC,  and  Mr.  Charles  C.  Calhoun,  Jr., 
Assistant  Director,  CERC;  and  under  the  direct  supervision  of  Mr.  C.  E, 
Chatham,  Chief,  WDD,  and  Mr.  Davidson,  Chief,  Wave  Research  Branch,  WDD.  The 
models  were  operated  by  Mr.  Robert  L.  Tingle,  Jr.,  Laboratory  Technician,  WDD. 

During  publication  of  this  report.  Dr.  Robert  W.  Whalin  was  Director  of 
WES.  COL  Leonard  G,  Hassell,  EN,  was  Commander. 


1 


CONTENTS 


Pap.e 

PREFACE  .  1 

CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC)  UNITS  OF  MEASUREMENT  .  3 

PART  I:  INTRODUCTION .  4 

Background  .  4 

Purpose .  .  .  4 

PART  II;  DEFINITION  OF  TEST  ?.\RAMETERS .  5 

PART  III:  THE  MODEL .  11 

Test  Facility .  11 

Test  Structure .  !  1 

Sounding  Equipment  .  16 

Signal  Generation  .  16 

PART  IV:  TEST  PROCEDURES  AND  TEST  RESULTS .  18 

Selection  of  Test  Conditions .  18 

Test  Procedures .  18 

Test  Results .  19 

PART  V:  DISCUSSION  AND  ANALYSIS .  21 

Runup . 21 

Damage  Level  .  27 

PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS .  42 

REFERENCES .  44 

APPENDIX  A:  DATA  TABLES . Al 

APPENDIX  B:  EXAMPLE  PROBLEMS . Bl 

APPENDIX  C:  NOTATION  .  Cl 


2 


CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement 
(metric)  units  as  follows; 

_ Multiply _ 

cubic  feet 

feet 

inches 

pounds  (mass) 

pounds  (mass)  per 
cubic  foot 

square  inches 


used  in  this  report  can  be 


_ By 

0.02831685 

0.3048 

2.54 

0.4535924 

16.01846 

6.4516 


converted  to  SI 

To  Obtain 

cubic  metres 

metres 

centimetres 

kilograms 

kilograms  per 
cubic  metre 

square  centimetres 
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USE  OF  A  RUBBLE  BERM  FOR  REDUCING  RUNUP.  OVERTOPPING. 
AND  DAMAUE  ON  A  IV  TO  2H  RIPRAP  SLOPE 


Experimental  Model  Investigation 

PART  I :  INTRODUCTION 

Background 

1.  Erosion  of  exposed  soil  embankments  by  waves  and  currents  is  a 
serious  problem  along  coastal  and  inland  shores.  A  protective  revetment  of 
graded  quarry  stone  or  "riprap"  is  commonly  used  to  provide  shore  protection 
because  of  its  relatively  low  cost,  durability,  and  availability,  and  because 
the  roughness  and  porosity  of  the  stone  is  effective  in  dissipating  wave 
energy  and  runup. 

2.  Rising  water  levels,  larger  boat  wakes,  and  increasing  values  of 
land  protected  by  revetments  may  necessitate  improvements  to  a  revetment  to 
increase  the  protection  provided.  One  method  of  improving  the  performance  of 
a  revetment  is  to  add  an  attached  berm  in  front  of  the  revetment.  Unfortu¬ 
nately,  there  is  very  little  design  guidance  on  the  use  of  revetments  with 
fronting  berms. 


Purpose 

3.  The  purpose  of  this  investigation  was  to  develop  design  guidance  on 
methods  of  reducing  runup  and  overtopping  of  revetments  based  on  data  col¬ 
lected  from  laboratory  tests  of  wave  runup  and  overtopping  on  riprap  revet¬ 
ments  with  a  slope  of  1:2  (IV  to  2H) .  Tests  were  conducted  in  a  wave  flume 
with  spectral  capabilities  at  the  Coastal  Engineering  Research  Center  (CERC) 
of  the  US  Army  Engineers  Waterways  Experiment  Station  (WES)  in  Vicksburg,  MS. 
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PART  II;  DEFINITION  OF  TEST  PARAMETERS 


4.  Inconsistencies  among  authors  in  notations,  definitions  of 
parameters,  and  the  methods  by  which  a  value  for  a  parameter  is  obtained 
greatly  complicate  the  task  of  comparing  results  from  different  studies.  In 
this  report,  notations  will  follow  guidelines  published  by  the  International 
Association  for  Hydraulic  Research  in  its  "List  of  Sea  State  Parameters" 
(1986).  Additional  parameters,  definitions,  and  method  used  to  determine  the 
value  of  certain  parameters,  are  given  below. 

5.  Wave  heights  used  in  this  report  are  the  heights  of  the  zeroth 

moment  (H^o)  obtained  as  four  times  the  square  root  of  the  zeroth 

moment  of  the  potential  energy  spectrum.  The  incident  spectra 

are  separated  from  the  H^o's  of  the  reflected  spectra  by  the  method  of  Goda 
and  Suzuki  (1976) ,  using  a  three-gage  array.  Two  arrays  are  used,  one  to 
measure  the  H^^'s  near  the  wave  board  (Array  1)  and  one  near  the  structure 
toe  (Array  2) . 

6.  Peak  period  (Tp^  is  the  wave  period  associated  with  the  highest 
energy  density  of  the  spectrum.  This  Tp  was  obtained  by  dividing  the 
spectrum  into  256  bands  and  taking  the  reciprocal  of  the  midpoint  frequency 
causing  the  highest  energy  density  over  11  adjacent  bandwidths . 

7 .  Wave  heights  and  periods  are  frequently  reported  in  other  investiga¬ 

tions  in  terms  of  significant  wave  height  (H^)  and  average  wave  period  (T^)  , 
where  is  the  average  of  the  one- third  highest  waves.  Both  and  T^ 

are  included  in  the  data  in  this  report  to  simpiity  comparison  to  other  inves¬ 
tigations.  Average  wave  periods  in  this  report  were  determined  as 


T 


z 


(1) 


where  m^,  and  m2  are  the  zeroth  and  second  moments  of  the  potential  energy 
spectrum,  respectively. 

8.  The  spectral  width  or  peakedness  determined  from  the  wave  record  is 
given  as  Qp  ,  defined  by  Goda  (1970)  as 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  C) . 
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0.  = 


f  flS(f)]'^  df 


<2) 


where  f  is  frequency  and  S(f)  is  the  wave  spectral  density  function  for 
the  given  frequency. 

9.  The  surf  parameter,  the  ratio  of  structure  slope  to  square  root  of 
wave  steepness,  is  frequently  used  as  an  indicator  of  wave  conditions  at  the 
structure  (Battjes  1974)  and  as  a  means  of  nondimens ionalizing  the  wave 
period.  The  surf  parameter  is  defined  as 


tan  a 


2  n  H 

gl^ 


1/2 


(3) 


or,  equivalently. 


5  = 


tan  a 


(4) 


where 

^  -  surf  parameter 

tan  a  ->  tangent  of  the  angle  of  the  revetment  to  the  horizontal 
g  =  acceleration  of  gravity 

Lg  -  deepwater  wavelength  determined  from  the  wave  period,  T 
The  wave  height  (H)  and  period  used  to  determine  ^  in  this  report  are 
and  Tp  .  The  physical  rationale  for  this  parameter  is  discussed  in  Battjes 
(1974) . 

10.  The  reflection  coefficient  is  commonly  defined  ..»s  the  ratio  of 
reflected  wave  height  to  incident  wave  height.  This  definition  is  clearly 
inappropriate  when  reflected  and  incident  wave  neights  are  described  by 
different  spectra.  Reflection  coefficients  were  therefore  determined  by  the 
energy  of  the  respective  spectra,  following  the  method  of  Goda  and  Suzuki 
(1976). 


(5) 
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where  Kj.  is  the  reflection  coefficient  and  Ep,  and  Ej  are  the  energy  of 
the  reflected  and  incident  spectra,  respectively. 

11.  Reflected  wav,  aeight  is  obtained  as  the  product  of  the  reflection 
coefficient  and  the  i'-ident  wave  height. 

H,  =  K,  (6) 

where  is  the  reflected  wave  height. 

12.  Runup  (R^ax)  defined  as  the  vertical  distance  above  the  still- 
water  level  (swl)  that  a  wave  surges  up  the  revetment,  or  the  upper  limit  of 
"green"  water  since  it  does  not  include  spash  or  spray.  The  elevation  of 
maximum  runup  was  determined  visually  and  then  measured  with  a  point  gage. 

13.  The  berm  in  front  of  the  revetment  is  described  by  its  height  and 
width.  Berm  width  (Wg)  is  defined  as  the  distance  along  the  horizontal  top  of 
the  berm.  Height  of  the  berm  (hg)  is  defined  as  the  vertical  distance  from 
the  toe  of  the  revetment  without  the  berm  to  the  top  of  the  berm. 

14.  Figure  1  illustrates  a  typical  damage  profile  on  a  riprap  slope. 

An  area  of  erosion  (Aj)  is  seen  near  the  still-water  level,  with  the  stones 
displaced  from  the  area  of  erosion  being  deposited  on  the  lower  slope  (A3)  or, 
particularly  on  very  flat  slopes,  on  the  upper  slope  (Aj)  .  Damage  to  the 
revetment  was  determined  based  on  the  area  of  erosion. 

15.  To  obtain  the  respective  areas,  the  six  sounding  points  on  each  of 

the  horizontal  sounding  lines  (paragraph  23)  were  averaged  to  give  a  single 
cross-sectional  profile  of  the  slope.  The  points  on  the  profile  then  were 
connected  by  cubic  splines,  and  the  before-  and  after- testing  profiles 
compared.  The  area  of  erosion  was  determined  as  the  area  between  the  two 
curves  where  the  after- testing  curve  lay  below  the  before  -  testing  curve.  The 
area  was  calculated  by  integrating  between  the  two  lines.  The  damage  profile 
from  Tp  -  2.25-sec,  -  0.50  ft  wave  conditions  is  shown  in  Figure  2. 

16.  Damage  to  the  revetment  was  defined  by  two  methods:  maximum  per¬ 
pendicular  penetration  of  the  erosion  into  the  armor  layer  (e,,,,,^)  and  $2 
damage . 

S2  =  - (7) 


7 


BEFORE  TESTING 
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Lgure  1.  Typical  damage  profile  on  a  riprap  revetment 


HORIZONTAL  DISTANCE  (fl.) 

filter  layer  -  -  ARMOR  LAYER  BEFORE  TESTING  ARMOR  LAYER  AFTER  TESTING 


where  (0^)50  is  the  nominal  diameter  of  the  median  armor  scone  size,  defined 
as 


where  Ujg  is  the  median  armor  stone  weight  and  w^.  is  the  unit  weight  of 
the  armor  stone. 
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PART  III:  THE  MODEL 


Test  Facility 


L7.  All  tests  were  conducted  in  a  3.0-ft-wide*  by  150-fL-long  by  3.0- 
ft-deep  wave  flume  (Figure  3).  A  1:20  slope  was  installed  in  the  bottom  of 
the  flume  starting  36.5  ft  from  the  wave  board  and  extending  for  10  ft, 
followed  by  a  1:100  slope  extending  to  the  test  structure. 

18.  The  flume  was  divided  lengthwise  into  two  1.5-ft-wide  channels 
starting  100  ft  from  the  wave  board  and  extending  past  the  structure.  A  wave 
absorber  was  placed  in  one  channel  while  the  structure  was  placed  in  the  otner 
channel  (113  ft  from  the  wave  board)  thus  minimizing  the  amount  of  reflection 
in  the  flume.  An  array  of  three  wave  gages  was  centered  21.5  ft  in  front  of 
the  wave  board  to  monitor  the  generated  signal.  A  similar  array  was  placed  in 
the  side  of  the  flume  with  the  test  structure,  centered  lOA . 5  ft  from  the  wave 
board,  to  be  used  in  separating  the  incident  and  reflected  wave  trains.  Wire 
resistance  staff  gages  were  read  at  10  Hz  to  monitor  the  water  surface 
elevation. 

19.  The  wave  flume  was  equipped  with  a  piston- type  wave  generator 
powered  by  an  electro-hydraulic  pump  and  controlled  by  a  computer  -  generated 
signal . 


Test  Structure 


20.  The  test  structure  modeled  a  1:2  slope  of  an  impervious  substratum 
protected  by  a  filter  layer  and  a  layer  of  riprap  (Figure  4).  Sand  was  glued 
to  a  plywood  board  to  provide  the  necessary  roughness,  and  the  board  installed 
in  the  flume  to  represent  the  existing  slope.  A  0 . 07 - f t - thick  layer  of 
crushed  stone  averaging  0.04  oz/stone  was  used  for  the  filter  layer.  The 
armor  layer  was  0 . 26  -  ft- thick  and  constructed  of  a  crushed  limestone  with  a 
specific  gravity  of  2.70,  blocky  to  angular  shape,  and  gradation  of 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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Figure  3.  Plan  and  profile  views  of  wave  flume 
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Figure  4.  Test  structure  cross  section  including  both  berms  tested 


=  1.79 


(9) 


where  Dsj  and  0^5  are  the  85-  and  15-percentile  diameters,  respectively,  on 
a  grain  size  distribution  curve,  with  all  armor  stones  falling  within  the 


range 


<  W  <4 


where  W  is  the  weight  of  an  individual  stone.  Gradation  of  the  armor  stone 
used  in  this  test  series  is  shown  in  Figure  5. 

21.  Riprap  is  commonly  sized  by  either  the  W50  or  the  nominal 
diameter,  (Dn)50  •  Tli®  armor  layer  used  in  these  tests  had  a  Wjg  of  0.22  lb 
and  a  (0^)50  of  0.11  ft.  These  dimensions  correspond  to  a  design  section 
based  on  Hudson's  equation  (Hudson  and  Jackson  1962)  for  a  design  wave  height 
of  0.30  ft.  Hudson's  equation  is  given  as 


^50  = 


(5^-1)  cote 


where 

H  -  monochromatic  wave  height 
Krjj  ->  stability  coefficient 

Sj  “  relative  specific  gravity  defined  as  specific  gravity  of  armor 
stone  divided  by  specific  gravity  of  water 

cot  -  is  the  cotangent  of  the  revetment  slope 
The  stability  coefficient  for  graded  angular  quarrystone  for  breaking  waves  is 
taken  from  the  Shore  Protection  Manual  (SPM)  (1984)  as  2.2.  For  irregular 
wave  tests,  design  is  based  on  the  average  of  the  highest  10  percent  of  waves 
(Hjo) .  The  design  H^g  of  0.3  ft  corresponds  to  a  H,,^  of  about  0.24  ft. 

22.  The  filter  and  armor  stone  layers  were  placed  by  dumping  from  a 
small  shovel  in  a  manner  to  simulate  prototype  construction  with  a  bucket 
loader.  The  toe  of  the  structure  was  113  ft  from  the  wave  board. 
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Figure  5.  Armor  stone  gradation 


23.  After  placement,  the  filter  and  armor  layers  were  surveyed  by 
sounding.  The  sounding  apparatus  consisted  of  a  row  of  six  vertical  rods 
spaced  across  the  half  of  the  flume  width  containing  the  test  structure,  with 
the  outer  rods  located  0.36  ft  from  the  edge  of  the  flume  section  and  the 
center  rods  spaced  at  0.15- ft  intervals.  Each  of  the  rods  had  a  flat  circular 
foot,  0.10  ft  in  diameter,  attached  to  the  rod  by  a  ball  and  .socket  joint. 

The  rods  were  used  to  take  a  row  of  soundings  across  the  test  structure.  A 
row  of  soundings  was  taken  at  horizontal  intervals  of  0.10  ft  down  the  entire 
slope  of  the  revetment.  Soundings  were  taken  before  and  after  each  test. 


Signal  Generation 


24.  Software  was  developed  at  CERC  to  generate  a  spectral  signal  to 
control  the  wave  board.  A  signal  was  generated  to  simulate  a  JONSWAP  (Joint 
North  Sea  WAve  Project)  spectrum  for  each  peak  period  using  high-  and  low- 
frequency  cutoffs  of  3.0  percent  of  the  energy  spectrum,  high-side  decay 
factor  of  0.09,  low-side  decay  factor  of  0.07,  peak  enhancement  factor  of  3.3, 
and  a  wave  height  (H,ao)  least  as  large  as  the  largest  in  the  test 

series.  Different  H,no*s  for  the  tests  were  obtained  by  varying  the  gain  to 
the  electro-hydraulic  pump  controlling  the  wave  board.  Computer  limitations 
dictated  that  30  min  was  the  maximum  signal  length  that  could  be  stored; 
therefore  a  separate  signal  was  generated  for  each  30  min  of  a  test  run,  using 
a  different  random  seed  for  each  signal.  The  design  spectrum  is  compared  with 
the  generated  spectrum  in  Figure  6  for  the  first  30  min  period  of  the  Tp  - 
2.25-sec  test  series. 
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WINDOW  MODE  «  1  (COSINE  SQUARED)  WINDOW  FRACTION  -  0.10 

DETREND  MODE  -  1  (MEAN  REMOVED)  BANDS  AVERAGED  •=  10 

WAtfcR  DETREND  FUNCTION  -  (  0.2005E-03)  +  (  O.OOOOE+00)*T  +  (  0.0O0OE+00)*CT**2) 
BOARD  DETREND  FUNCTION  -  (  0,2325E-03)  +  (  O.OOOOE+00)*T  +  (  O.0OOOE+O0)*CT**2) 


UNIDIRECTIONAL  WAVE  TANK  PISTON-TYPE  GENERATOR  DIGITAL  DRIVER  RLE  HEADER 


OUTPUT  RLE  -  0T2T225H6.DAT 
GENERATOR  TYPE  -  PISTON 
TRANSFER  FUNCTION 
SPECTRUM  -  JONSWAP  SPECTRUM 
2.250  (SEC)  DESIGN  PEAK  PERIOD 
6.000  (INCHES)  DESIGN  HMD 
3.30  GAMMA 
0.0700  SIGMA-LOW 
0.0900  SIGMA-HIGH 

0.00310  COMPUTED  ALPHA 

7.489  INCHES  MAX.  POS.  BOARD  STROKE 
-8.228  INCHES  MAX.  NEG.  BOARD  STROKE 
0.3544  (HZ)  LOW  FREQ.  CUTOFF 
1.0096  (HZ)  HIGH  FREQ.  CUTOFF 


JOHN  AHRENS  CORRECTION  -  NO 
24.000  (INCHES)  WATER  DEPTH  AT  BOARD 
0.030  X100PCT  MO  LOW  FREQ.  CUTOFF 
0.970  X^OOPCT  MO  HIGH  FREQ.  CUTOFF 
1800  (SEC)  TIME  SERIES  LENGTH 
10  (STEPS/SEC)  UPDATE  RATE 
256  UNES  IN  SOURCE  SPECTRUM 

1989  -  RAND.  NO.  GEN.  SEED 
16.64  -  (IN**2/HZ)  DES.  H20  SMAX 

0.4444  -  (HZ)  DES.  PEAK  FREQ. 


Figure  6.  Design  and  simulated  spectra  for  first 
30  min  of  2.25-sec  wave  period  tests 
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PART  IV:  TEST  PROCEDURES  AND  TEST  RESULTS 


Selection  of  Test  Conditions 


25.  Tests  were  conducted  for  water  depths  at  the  toe  of  the  revetment 
of  0.66,  0.76,  and  0.86  ft,  with  offshore  water  depths  of  2.00,  2.10,  and 
2.20  ft,  respectively.  Peak  wave  periods  tested  included  1.25-,  2.25-,  3.0-, 
and  3.5-sec  waves,  and  wave  heights  of  the  zeroth  moment  for  each  period 
included  0.10,  0.20,  0.30,  0.40,  0.50,  0.60,  and  0.66  ft. 

26.  The  length  of  each  test  run  was  3,500  times  the  peak  period  of  the 
test,  or  until  failure  occurred,  whichever  came  first.  Failure  was  defined  as 
having  the  filter  layer  exposed.  If  failure  was  observed,  tests  with  that 
period  were  suspended  and  larger  H^'s  were  not  tested.  All  tests  were  con¬ 
ducted  with  a  JONSWAP  spectrum  (Hasselman  et  al.  1973). 

Test  Procedures 


27.  Each  test  was  run  in  cycles  of  approximately  30  min  each.  For 
example,  the  tests  with  a  2.25-sec  peak  wave  period  were  run  for  7,875  sec 
(3,500  times  2.25  sec),  or  131.25  min,  in  four  cycles  of  about  30  min  each 
plus  one  short  cycle  to  obtain  the  total  time  required.  Wave  data  were 
collected  for  approximately  500  waves  near  the  beginning  and  end  of  each  test. 
A  test  was  terminated  if  failure  occurred. 

28.  Wave  runup  was  observed  for  256  sec  near  the  middle  of  a  test  run. 
The  elevation  of  the  maximum  runup  observed  then  was  measured  with  reference 
to  a  specified  datum  by  use  of  a  point  gage. 

29.  Soundings  were  taken  before  and  after  each  test,  and  the  structure 
was  completely  rebuilt  whenever  significant  damage  was  observed.  The  pro¬ 
cedure  used  for  the  tests  is  given  below. 

a.  Place  and  survey  the  filter  stone. 

b.  Place  and  survey  the  armor  stone  (including  the  berm,  if 
used) . 

c.  Generate  waves  for  3,500  times  the  design  peak  period. 

(1)  Collect  wave  data  for  500  waves  near  the  start  of  the 
test. 
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(2)  Measure  the  runup. 

(3)  Collect  wave  data  for  500  waves  near  the  end  of  the  test. 
If  failure  is  expected,  collect  data  for  500  waves  prior 
to  the  expected  failure. 

d.  Survey  the  riprap  layer  and  berm. 

e.  Measure  the  water  in  the  overtopping  catch  basin. 

f.  Rebuild  the  structure. 

30.  If  the  damage  was  slight,  only  the  riprap  layer  and  berm  were  re¬ 
built,  and  the  procedure  for  testing  the  next  wave  height  was  begun  at  step  b, 
above.  The  structure  was  always  torn  down  and  rebuilt  from  the  pl3WOod  board 
before  the  next  series  of  tests  with  a  new  wave  period  was  conducted. 

Test  Results 


31.  Test  results  are  given  in  Appendix  A  as  Tables  Al,  A2,  A3,  A4,  and 
A5.  Tables  Al  and  A2  list  data  from  the  tests  of  the  revetment  without  berm; 
Tables  A3  and  A4  list  the  results  from  the  revetment  with  berm,  and  Table  A5 
lists  damage  to  the  structure  as  determined  by  the  soundings. 

32.  As  discussed  in  paragraph  27,  multiple  cycles  were  required  for 
each  test  run.  Each  cycle  was  labeled  consecutively  a,  b,  c,  etc.,  and  data 
were  collected  during  the  first  and  last  cycle.  For  example,  tests  with  a 
2. 25 -sec  peak  wave  period,  which  required  four  cycles,  had  data  collected 
during  cycles  a  and  d.  The  letter  following  the  run  number  in  Tables  Al 
through  A4  refers  to  the  cycle  during  which  the  data  were  collected.  The 
structure  was  surveyed  only  after  the  end  of  the  test  run.  Therefore  the 
damage  information  in  Table  A5  refers  only  to  the  run  number  and  does  not 
include  the  cycles.  Because  damage  to  the  berm  was  not  considered  as  damage 
to  the  revetment,  only  damage  values  for  the  upper  part  of  the  revetment 
(above  the  level  of  the  berm)  are  used  in  this  report  and  are  referred  to 
simply  as  $2  . 

33.  Scale  effects  in  hydraulic  model  investigations  are  an  important 
consideration,  but  one  that  is  difficult  to  address  due  to  conflicting  results 
from  different  studies.  For  the  range  of  tests  presented  here,  Hudson  and 
Davidson  (1975)  indicate  that  stability  numbers  in  the  model  tests  will  be 
conservative  by  about  10  percent  compared  with  prototype  results,  Burcharth 
and  Frigaard  (1988)  report  that  scale  effects  will  be  negligible,  and 
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Broderick  and  Ahrens  (1982)  show  that  scale  effects  in  the  filter  layer  may 
cause  the  runup  values  in  the  model  tests  to  be  greater  than  would  be  found  at 
prototype  scales.  Although  scale  effects  are  probably  present  in  this  data 
set,  it  seems  reasonable  to  assume  that  they  are  either  negligible  or 
conservative . 
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PART  V:  DISCUSSION  AND  ANALYSIS 


34.  Wave  runup  and  overtopping  data  were  collected  from  a  typical 
revetment  design  subjected  to  spectral  wave  attack.  As  a  means  of  reducing 
the  runup  and  overtopping  of  the  revetment,  a  berm  has  been  added  in  front  of 
the  revetment.  Comparisons  of  results  of  the  tests  with  and  without  the  berm 
will  be  used  to  develop  design  guidance  for  use  by  the  field  to  reduce  runup 
and  overtopping  of  rubble  revetments.  Structural  stability  and  reflection 
characteristics  also  are  being  documented. 


Wave  Runup 


35.  Ahrens  and  McCartney  (1975)  documented  the  influence  of  surf 
conditions  on  runup,  and  an  empirical  model  using  the  surf  parameter  to 
predict  maximum  runup  has  been  developed  by  Ahrens  and  Heimbaugh  (1988).  The 
empirical  model  is  given  by 

^  ^  (12) 

1 . 0  +  b  5 


where  is  maximum  wave  runup  and  a  and  b  are  dimensionless  runup 

coefficients  determined  from  regression  analysis.  Values  of  a  and  b  were 
determined  as  1.022  and  0.247,  respectively. 

36.  Analysis  of  the  data  collected  in  this  study  indicates  that  two 
additional  terms  are  needed  to  account  for  the  influence  of  a  berm.  These 
additional  terms  account  for  the  height  and  width  of  the  berm.  Through  data 
analysis,  it  was  found  that  the  most  effective  variables  were  relative  berm 
width,  Wb/(H,„„L„)''^  ,  and  relative  berm  height,  hg/d,  .  Multivariable  analysis 
was  used  to  determine  the  following  equation  to  predict  the  relative  runup, 


Rmax 

H„o 


exp 


r  +  r  1  ' 

+  r 

r  M  ^  J 

^2 

1  1/2 

(13) 


where  Co  ,  C,  ,  and  Cj  are  dimensionless  regression  coefficients  given  by 
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Co  -  0.669 
Cl  =  -10.4 
C2  -  -0.152 

37.  Equation  13  shows  liL._le  or  no  systematic  error  in  predicting 

as  shown  by  error  analysis  in  Figures  7,  8,  and  9.  Figure  7  gives  the  ratio 
of  the  predicted  to  the  observed  as  a  function  of  the  relative 

berm  width.  Figure  8  has  the  same  ordinate  as  Figure  7  plotted  versus  the 
relative  berm  height.  Figure  9  is  an  enlargement  of  the  portion  of  Figure  8 
for  tests  with  berms  so  that  the  data  point  symbols  can  be  more  readily 
identified.  In  Figures  7,  8,  and  9,  all  data  points  are  denoted  by  an  integer 
based  on  the  relative  wave  size,  H^/(Dn)5o  .  Since  deterioration  of  the  berm 
increases  with  increasing  wave  size  relative  to  the  stone  size,  the  purpose  of 
using  the  relative  wave  size  in  the  figures  is  to  see  if  berm  derioration 
affects  the  prediction  of  .  Figures  7,  8,  and  9  indicate  there  is 

little  systematic  error  in  predicting  R^^^  for  tests  with  a  berm. 

38.  Figure  10  shows  relative  runup  versus  wave  steepness,  ,  for 

tests  without  a  berm.  On  a  planar  riprap  slope  without  a  berm,  wave  steepness 
is  the  most  important  variable  influencing  relative  runup.  Also  shown  in 
Figure  10  is  the  portion  of  Equation  13  without  the  berm  effects  correction 
factor ,  i . e . , 


^^Toax 

Hn,o 


exp 


Co  ^  C, 

(  H  \1 

*^mo 

L 

V  •^0  n 

(14) 


where  Cq  and  are  the  same  as  given  with  Equation  13.  For  tests  without 

berms,  the  berm  effects  term  in  Equation  13  drops  out  to  leave  Equation  14. 

In  Figure  10,  it  can  be  seen  that  the  data  scatter  around  Equation  14  appears 
random,  indicating  a  lack  of  systemic  error  associated  with  wave  steepness. 
Figure  10  also  shows  the  runup  equation  developed  by  Ahrens  and  Heimbaugh 
(1988)  for  plane  riprap  slopes  (Equation  12).  In  Figure  10  it  can  be  seen 
that  Equation  12  overpredicts  the  relative  runup  but  follows  the  data  trend 
well.  It  is  believed  that  this  difference  is  caused  by  systematic  differences 
in  the  visual  readings  of  R^^,,  between  ob*^«rvers, 

39.  The  portion  of  Equation  13  containing  the  berm  characteristics  can 
be  thought  of  as  a  runup  reduction  factor  that  can  be  defined  by 
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Figure  7.  Ratio  of  predicted  to  observed  runup  as  a  function  of  relative  bertn  width. 

Integers  denote  wave  size,  H„„/(D„)5g 


Figure  8.  Ratio  of  predicted  to  observed  runup  as  a  function  of  relative  berm 

Integers  denote  wave  size,  H^/(Dn)5o 
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Figure  9.  Section  taken  from  Figure  8  enlarged  to  show  detail 


r 


=  exdc 


W. 


♦ 


(15) 


Figure  11  shows  lines  of  constant  r  as  a  function  of  the  relative  berm  width 
and  relative  berm  height  over  the  range  of  data  collected  in  this  study  based 
on  the  value  of  Cj  -  -0.152  .  Figure  11  helps  to  visualize  the  joint 
influence  of  berm  width  and  berm  height  in  reducing  wave  runup  and  also 
demonstrates  that  Equation  15  has  a  logical  functional  form. 

40.  To  further  generalize  the  findings  from  this  study,  Equations  12 
and  13  were  combined  in  the  following  manner 

(16) 

H„«  1  *  hi 


The  "a"  runup  coefficient  was  dropped  from  Equation  12  since  it  is  very 
close  to  1.0  and  replaced  with  r  ,  which  has  a  limiting  value  of  1,0  for 
plane  slope  revetments.  Figure  12  is  a  scatter  plot  of  predicted  using 

Equation  16  versus  observed  .  It  is  clear  from  Figure  12  that 

Equation  16  follows  the  trend  of  the  data  well  but  provides  conservative 
estimates  of  R^  .  As  noted  in  paragraph  38,  systemic  differences  between 
runup  observations  in  this  study  and  those  in  Ahrens  and  Heimbaugh  (1988)  may 
account  for  the  differences  in  observed  and  predicted  runup  in  Figure  12. 


Damage  Level 


41.  In  addition  to  reducing  wave  runup  and  overtopping,  the  rubble  berm 
also  reduced  damage  levels  on  the  revetment  slope  above  the  berm.  Figures  13a 
through  13d  show  the  progressive  deterioration  of  the  berm  with  increasing 
wave  heights.  The  protection  furnished  to  the  revetment  is  obvious,  as  seen 
when  these  figures  are  compared  with  Figures  14a  through  14d.  Figures  14a 
through  14d  show  the  progressive  deterioration  of  a  revetment  without  a  berm 
after  being  subjected  to  the  same  wave  conditions  used  in  Figures  13a  through 
13d. 

42.  Quantification  of  the  improvement  can  be  presented  in  a  manner 
similar  to  the  method  used  for  wave  runup  in  Figure  11,  i.e.,  by  using  a 
reduction  factor.  Damage  to  the  structure  above  the  berm  can  be  quantified 
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Figure  11.  Runup  reduction  factor,  r  ,  as  a  function  of  relative  berm  height 
and  width,  using  Cj  -  -0.152  in  Equation  15 
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Maximum  runup  predicted  by  Equation  16  versus  observed  maximum  runup 
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either  by  defining  the  damage  for  the  region  above  the  berm,  or  by  the 
depth  of  erosion  into  the  armor  layer,  e,,,,  ,  The  following  equations  were 
developed  using  multivariable  analysis  to  predict  damage  levels. 

Sj  =  1.0  exp 


W, 


(H„oLo) 


n 


(17) 


=  Cc 


exd 


w. 


(18) 


where  N,  is  the  stability  number,  defined  below,  t,  is  the  armor  layer 
thickness,  and  Cj  through  C,  are  dimensionless  regression  coefficients 
given  by 

Cj  -  2.96 
C4  -  -3.56 
Cj  -  0.104 
C*  -  2.08 
C-,  -  -2.73 

The  stability  number  is  defined  as 


Ns  = 


w 


1/3 


H„.c 


W, 


1/3 

50 


(Sj  -  1) 


(19) 


43.  Equations  17  and  18  were  developed  to  show  how  the  reduction  factor 
concept  introduced  for  runup  can  also  be  applied  to  damage  levels  and  are  not 
necessarily  recommended  for  prediction  purposes  since  data  scatter  is  very 
high  in  the  damage  variables.  However,  all  terms  in  the  equations  are  highly 
significant  and  the  equations  seem  to  reflect  the  trends  in  the  data. 

44 .  Damage  can  be  defined  as  follows ; 


r.  =  exdC 


Wp 


I  (H^Lo) 


1/3 


(20) 


=  exp^ 


Wp 


(21) 
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where  and  are  the  reduction  factors  for  S2  damage  and  e^aj; 

damage,  respectively.  Figures  15  and  16  show  curves  of  constant  reduction  as 
functions  of  relative  berm  width  and  relative  berm  height  based  on 
Equations  20  and  21.  These  figures  help  visualize  the  influence  of  a  berm  in 
reducing  damage  to  a  revetment. 
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Figure  15.  S2  damage  reduction  factor,  r,  ,  as  a  function  of  berm  height  and  width 
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Relative  Berm  Width,  WB/(Hmo*Lo)’« 


PART  VI;  CONCLUSIONS  AND  RECOMMENDATIONS 


45.  This  study  developed  an  equation  (Equation  16)  to  predict  the  upper 
limit  of  irregular  wave  runup  on  planar  riprap  revetments  or  on  riprap  revet¬ 
ments  fronted  by  a  rubble  berm.  The  equation  considers  both  the  height  and 
width  of  the  berm  and  provides  a  way  to  estimate  the  effectiveness  of  a  berm 
in  reducing  wave  runup.  Although  the  berm  was  only  tested  with  a  1:2  slope, 
the  reduction  factor  is  believed  to  be  applicable  to  other  slopes  when  used  in 
Equation  16 .  Since  Equation  16  can  be  used  to  calculate  the  runup  on  both 
bermed  and  plane  riprap  revetments,  it  is  superior  to  the  method  proposed  in 
Repair,  Evaluation,  Maintenance,  and  Rehabilitation  Research  (REMR)  Program 
Technical  Note  CO-RR-1.3,  Supplement  3,  which  provided  only  a  reduction  factor 
for  runup  due  to  the  influence  of  a  berm.  In  addition,  Equation  16  is  based 
on  a  larger  data  set  and  subsequent  analysis  beyond  that  given  in  REMR 
Technical  Note  CO-RR-1.3.  An  updated  version  of  REMR  Technical  Note  CO-ElR-1.3 
will  be  available  in  Supplement  6  to  reflect  these  changes. 

46.  A  rubble  berm  not  only  reduces  wave  runup  but  also  increases  the 

stability  of  the  revetment.  The  berm's  influence  on  stability  is  considerably 
greater,  as  a  percent,  than  the  reduction  in  runup;  i.e.  a  modest  reduction  in 
runup  corresponds  to  a  substantial  improvement  in  stability.  Typical  reduc¬ 
tions  in  runup  due  to  a  berm  observed  during  this  study  were  in  the  range  of  5 

to  15  percent.  It  is  also  believed  that  a  modest  reduction  in  runup  will 

translate  into  a  large  reduction  in  overtopping  rates.  These  intuitive  con¬ 
cepts  can  be  supported  by  making  a  few  logical  assumptions  and  calculating  the 
consequences . 

47.  The  reduction  in  overtopping  can  be  estimated  by  using  the  method 

shown  in  the  SPM  and  discussed  further  by  Weggel  (1976)  .  Assume  that  R,„„  - 
1.15  *  F;  i.e.  the  potential  maximum  runup  is  15  percent  greater  than  the 
freeboard  of  a  plane  slope  riprap  revetment.  To  reduce  the  overtopping  rate, 
a  berm  is  added  to  the  revetment  that  reduces  the  maximum  runup  10  percent  or 

R,„„  -  (1.15  *  0.90)  *  F  -  1.035  *  F,  where  F  is  the  structure  freeboard. 

This  rather  modest  reduction  in  runup  causes  a  reduction  in  the  overtopping 
volume  by  a  factor  of  about  8.8  if  a  value  of  Weggel 's  overtopping  parameter 
alpha  -  0.07  is  used  (see  Example  1,  Appendix  B) .  These  reductions  are  in  the 
volume  of  water  per  unit  length  of  structure  for  the  maximum  potential  runup . 
The  average  overtopping  rates  could  be  expected  to  be  reduced  more  than  the 
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maximum  rates  since  the  berm  would  allow  fewer  waves  to  overtop.  Values  of 
alpha  selected  are  reasonable  for  a  riprap  revetment  based  on  estimating  over¬ 
topping  rates  using  the  potential  runup  approach  given  by  Weggel  and  the  SPM. 

48.  Reductions  in  damage  for  modest  reductions  in  runup  are  equally 
impressive  as  the  reductions  in  overtopping  rates.  Compare  the  reduction 
factors  for  runup  and  Sj  damage  gi\an  by  Equations  15  and  19  for  a  berm  pro¬ 
ducing  a  10-percent  reduction  in  runup.  The  damage  reduction  factor  cor¬ 
responding  to  r  -  0.90  is  r,  -  0.085  (see  Example  2,  Appendix  B) .  This 
comparison  indicates  that  a  10-percent  reduction  in  runup  will  reduce  Sj 
damage  by  a  factor  of  almost  12. 

49.  Although  comparable  reductions  in  runup  and  overtopping  may  be 
obtained  with  a  smaller  volume  by  increasing  the  crest  height  of  the  revetment 
rather  than  adding  a  berm,  the  use  of  a  berm  should  prove  of  value  in  loca¬ 
tions  where  an  increased  crest  height  is  not  possible  or  undesirable.  The 
advantages  of  increased  protection  provided  by  the  berm  should  also  be 
considered. 

50.  Maximum  runup  elevations  from  this  study  are  based  on  visual 
observations.  These  observations  appear  to  be  consistently  lower  than  similar 
observations  made  during  other  studies  (Ahrens  and  Heimbaugh  1988) .  For  this 
reason,  data  analysis  was  biased  towards  a  conservative  interpretation  of  the 
findings  from  this  study.  Because  of  the  difficulty  in  achieving  consistent 
results  among  different  observers  of  maximum  runup  elevations,  it  is  recom¬ 
mended  that  an  electronic  runup  gage  be  developed.  A  runup  gage  would  not 
only  give  consistent  results  from  study  to  study,  but  could  also  provide  more 
information  about  irregular  wave  runup  elevations  on  rough  and  porous  slopes. 
With  a  gage,  it  is  anticipated  that  statistically  stabler  runup  parameters 
could  be  computed  that  would  yield  better  insight  into  the  ability  of  a  berm 
to  improve  the  performance  of  a  rubble  structure.  Using  a  maximum  value 
contributed  to  the  considerable  data  scatter  in  this  study  and  the  difficulty 
in  interpreting  results. 
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Table  5 

Damage  Sustained  by  the  Revetment  During  Test  Runs 


w  # 

UPPER 

S2 

DAMAGE 

MAXIMUM 

PERPEN¬ 

DICULAR 

ARMOR 

EROSION 

1 

0.560 

0.010 

2 

0.883 

0.024 

3 

0.558 

0.013 

4 

1.323 

0.053 

5 

2.161 

0.049 

6 

6.517 

0.112 

7 

11.473 

0.121 

8 

0.514 

0.008 

9 

1.252 

0.036 

10 

5.538 

0.090 

11 

11.219 

0.157 

12 

21.797 

0.215 

17 

0.458 

0.006 

18 

0.458 

0.006 

19 

0.770 

0.036 

20 

1.220 

0.030 

21 

3.637 

0.085 

22 

9.358 

0.126 

23 

8.765 

0.158 

24 

0.338 

0.007 

25 

1.523 

0.037 

26 

12.590 

0.137 

27 

25.188 

0.191 

28 

0.486 

0.017 

29 

2.442 

0.044 

30 

12.220 

0.171 

32 

3.258 

0.097 

33 

12.191 

0.146 

34 

17.883 

0.195 

35 

1.240 

0.030 

36 

37.603 

0.286 

UPPER 

MAXIMUM 

PERPEN¬ 

DICULAR 

RUN  # 

S2 

ARMOR 

DAMAGE 

EROSION 

37 

0.189 

0.006 

38 

16.858 

0.194 

40 

0.257 

0.008 

42 

0.495 

0.014 

44 

0.927 

0.018 

47 

0.244 

0.008 

48 

4.881 

0.076 

49 

2.479 

0.048 

50 

0.301 

0.004 

51 

6.491 

0.088 

52 

1.758 

0.058 

53 

0.283 

0.006 

55 

0.209 

0.004 

57 

0.086 

0.004 

59 

1.946 

0.036 

60 

6.944 

0.129 

61 

9.695 

0.125 

63 

0.624 

0.010 

64 

0.693 

0.015 

65 

1.510 

0.041 

67 

0.973 

0.036 

68 

1.306 

0.054 

69 

6.841 

0.123 

70 

9.876 

0.104 

71 

17.981 

0.175 

85 

2.152 

0.044 

8  6 

1.171 

0.030 

88 

0.890 

0.021 

89 

0.730 

0.017 

90 

0.538 

0.016 

91 

1.610 

0.047 
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APPENDIX  B:  EXAMPLE  PROBLEMS 


B1 


EXAMPLE  1 


GIVEN:  Potential  runup  on  a  revetment  is  15  percent  higher  than  the 

freeboard. 

FIND:  Reduction  in  overtopping  rate  if  a  berm  reduces  the  runup  by  10  percent. 

SOLUTION:  The  current  method  given  in  the  Shore  Protection  Manual  (SPM) 

(1984)*  to  calculate  overtopping  rates  is 


*  3 

^6^0  «o 


11/2 


0.1085  . 
exp  -  In 


R  +  h  -  d 


R  -  h  +  d 


(Bl) 


where 


Q 

and  a 


h; 


R 

h 


overtopping  rate  per  unit  structure  length 

empirically  determined  coefficients  based  on  incident  wave 
conditions  and  structure  geometry 

equivalent  deep  water  wave  height 

maximum  potential  runup 

height  of  structure  crest  along  the  bottom 
depth  at  the  structure  toe 


This  may  be  rewritten  as  (Ahrens  and  Heimbaugh  1986) 


R  -  F 

1  0.1085 

l<8  <’5 J  1 

R  +  F 

J  « 

(B2) 


where  F  is  the  freeboard  of  the  structure. 

Given  that  the  potential  runup  without  a  berm  is  15  percent  greater  than 
the  freeboard,  then 

*^No  Berm  "  1-150F 


If  a  berm  reduces  the  runup  by  10  percent,  then 

RBerm  “  1-150F  *  0.90  -  1.035F 
The  effect  of  the  berm  on  the  overtopping  rate  is  given  by 


^No  Berm 
^Berm 


1.150F 

1.150F 


1.035F 

1.035F 


+ 


+ 


0.1085 

a 


p.  0.1085 

fJ  “ 


(B3) 


*  See  References  at  the  end  of  the  main  text. 


or, 


^Ko  Berm  _ 
^Betw 


I 0.150F\ 
\ 2.150F/ 
I 0.035F\ 
2.035F/ 


0.1085 


For  a  -  0.07,  the  berm  has  reduced  the  overtopping  rate  by  a  factor  of  8.8. 


EXAMPLE  2 


GIVEN; 


FIND: 

SOLUTION: 


A  given  berm  produces  a  10-percent  reduction  in  wave  runup,  i.e. 
r  -  0.90. 

What  is  the  damage  reduction  factor  (tg)? 


r  =  0.90  =  exdc 


(15  bis*) 


X.  =  ?  =  expcj— - 


(20  bis**) 


Since  Wg  ^  Hg  ,  Hj^q  ,  Lq  ,  and  dg  are  equal  and  C2  -  -0.152  (pp  25  of 
text)  and  -  -3.56  (pp  36  of  text),  the  above  reduces  to: 


r  ^  ^2 ' 

—  =  exd  — 

^3 


In  0.90 
In  (rj 


-0 ■ 152 
-3.56 


-0.105 
In  (rj 


-0.152 

-3.56 


*  See  p  27  in  the  main  text. 
**  See  p  38  in  the  main  text. 


B3 


-2.459 


In  r. 


APPENDIX  C :  NOTATION 


Cl 


kl 

A3 

alpha 

b 

Co 

Cl 


C6 

C7 

Di5 

^85 


'max 

% 

f 

F 

g 

H 

Hmo 

Hr 

Hs 

Kr 

Krr 

L, 


o 


m. 


o 


Regression  coefficient 

Area  of  accretion  above  the  still  water  level 
Area  of  erosion 

Area  of  accretion  below  the  still  water  level 

Overtopping  parameter 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

Regression  coefficient 

15-percentile  diameter  on  a  grain  size  distribution  curve 

85-percentile  diameter  on  a  grain  size  distribution  curve 

Nominal  diameter  of  the  median  stone  size 

Depth  at  structure  toe 

Energy  of  the  incident  wave  spectrum 

Maximum  perpendicular  penetration  of  erosion  into  armor  layer 

Energy  of  the  reflected  wave  spectrum 

Wave  frequency 

Structure  freeboard 

Gravitational  acceleration 

Wave  height 

Berm  height 

Wave  height  of  the  zeroth  moment 
Reflected  wave  height 

Average  wave  height  of  the  one-third  highest  waves 
Reflection  coefficient 
Stability  coefficient 
Deepwater  wavelength 

zeroth  moment  of  the  potential  energy  spectrum 


C2 


r 


s 

^ax 


S(f) 

S,. 


z 

W 


*'50 

Wn 


a 

7r 


Second  moment  of  the  potential  energy  spectrum 

Stability  number 

Spectral  width  or  peakedness 

Runup  reduction  factor 

Erosion  reduction  factor 

S2  damage  reduction  factor 

Maximum  vertical  height  above  still  water  level  of  wave  runup 

Damage  relative  to  size  of  armor  unit 

Spectral  density  function 

Relative  specific  gravity 

Wave  period 

Wave  period  associated  with  the  peak  energy  density 

Average  wave  period 

Weight  of  an  individual  armor  unit 
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